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Abstract 



O \ The formal Heisenberg equations of the Federbush model are linearized and then are 



directly integrated applying the method of dynamical mappings. The fundamental 
role of two-dimensional free massless pseudo-scalar fields is revealed for this procedure 
'sj" ' together with their locality condition taken into account. Thus the better insight into 

\ solvability of this model is obtained together with the additional phase factor for its 

. ' general solution, and the meaning of the Schwinger terms is elucidated. 

'NT 

> 
in 



1 Introduction 



Generalizing the methods, giving solutions in sufficiently simple models one can hope to 
develop a technique to be effective for the comparatively non-trivial and physically impor- 
<^ . tant problems [20] • Examples of such partly or entirely solvable models are models with 

four-fermion interactions in two-dimensional space-time and related to them non-linear 
bosonic models of the sine-Gordon (SG) type. The relentless interest in these models is 
due to the fact that their non-Abelian four-dimensional analogs are more or less success- 
r~| ! fully used for the analysis and explanation of various non-perturbative effects in modern 

theory of strong interactions, such as: description of the processes of quark hadronization 
and the phenomenon of spontaneous symmetry breaking |Hj, jSO]. They serve also as a 
"testing ground" for various non-perturbative methods [J, 0, p], [HI, [S]. As a rule an 
\ operator solution for these models is obtained with the help of proper operator ansatz 

[TT] . [T7] . [TH] . pn] instead of consecutive integration the corresponding equations of 
motion. 

Recently Faber, Ivanov in a series of papers (see e.g. [Uj, jlOp. following Morchio 
et al. ^H]) re-examined some ambiguities jZj of the Thirring model and elucidated the 
importance of existence of two-dimensional free massless scalar fields for its solution. 

One of the aims of the present work is to show a similar role of free massless pseudo- 
scalar fields for a direct step-by-step integration of the Heisenberg equations (HEqs) of the 
Federbush model and to advocate as a general method the corresponding linearization 
procedure, successfully used to solve some non-relativistic and relativistic phenomenologi- 
cal models ^Hl) HZ]- It is interesting to note that unlike free massless scalar field, the free 
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massless pseudo-scalar field has a well-defined generating functional JHl and appears in a 
definite sense as a gauge invariant object [2]. 

Another goal is to understand the role of Schwinger terms for this model. We show that, 
in spite of a long story of this sloppy question [221, ES], [2Hj the suggested linearization 
procedure allows naturally to adopt an idea of solution found in QED [21] . 

In order to linearize the HEqs we use the notion of physical fields [^, [201 • 
us recall that the HEqs are formal relations between operators, as long as they are not 
defined on a corresponding vector space. This implies, that in order to give a physical 
meaning to the quantum- field-system description by means of the Heisenberg fields ^{x), 
it is necessary to represent them in the space of physical states, that, in turn, necessitates 
to express them in terms of physical fields ip{x). The relationship between Heisenberg 
and physical fields is called dynamical mapping (DM), and may be formally expressed 
generally only as a weak equality: '^{x) = T[?/^(x)]. After quantization of the physical 
fields, satisfying some free field equations, one obtains the HEqs and DM in the normal- 
ordered form with respect to creation and annihilation operators of a chosen set of physical 
particles. This normal-ordered form in fact is equivalent in the sense of DM to their 
primary representations as a weak limit, but is, of course, representation dependent [1]. 

Since dynamical HEqs contain time derivatives it is necessary to choose initial con- 
ditions. There are two physically different choices leading respectively to two essen- 
tially different sets of physical fields: for asymptotic m-fields [H], [2015 for t — oo, 
one has lim ^{x^,t) = T['0j„(x^, — oo)]; alternatively for Schrodinger fields ^] one has 

t—r — OO 

lini^(a;^,t) = T['0(x^,O)] as t — > 0. The latter choice will be used here. It was advocated 

in our previous works [23 as some generalization of the usual interaction represen- 
tation [201 shown its usefulness on the bound-state problem, and in fact was already used 
for the construction of the solution of the Federbush model jllj . 

The work is organized as follows. First in section 2 we discuss a linearization procedure 
for HEqs. Then in section 3 we obtain their solutions in terms of Feynman functional 
integrals calculated in sections 4 and 5 by pure algebraic operator methods using various 
bosonization rules for the free fermion currents only and the properties of free massless 
pseudo-scalar fields only. It should be emphasized in order to avoid some confusions, that, 
as is shown below, the bosonization rules in fact could not be applied here to the local 
Heisenberg currents directly, and that in two-dimension space-time both the used in section 
4 bosonization rules: ()4.1|) for the local free massive fermion current and ()4.6|) for the local 
free massless fermion current, are immediate well-known consequence of conservation of 
both the vector currents ^H], [211. These two rules introduce however two different pseudo- 
scalar fields: the free massless field (/'(x) in (|4.6|) and the field ^(x) in (|4.H) which in fact is 
the sine-Gordon field, as is finally checked in subsection 4.2. These fields certainly serve as 
building blocks for the construction of solution, nevertheless, are necessary but insufficient 
ingredients to obtain the explicit operator solution of the model, as is shown in sections 4 
and 5. The role of Schwinger terms is also discussed in the Summary. 

2 Step 1. Linearization of Heisenberg equations 

The model under consideration describes a parity non-conserving local interaction of two 
different two-component fermion fields ^^(a;) in two-dimensional space-time with different 
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masses m^, ^ = ±1, and is defined by a Lagrangian density, bilinear in the two Noether 
vector currents J^(^q,^{x) = ^'^(x)7'^^'^(x). Splitting the total Hamiltonian of the model 
into free and interaction parts as 



oo 

€=±1 



(2.1) 
(2.2) 



and using the equal-time anticommutation relations for the Heisenberg operators^ 

{M/5,„(x),vI/^,,^(y)} =0, [^i,,^{x),^\,p{y)] = 5^^,6af, S{x^ - y^), (2.3) 



we get the formal HEqs for the fields ^^{x) 

ido^^ix) = [^^{x),H^^){t)] = [E^{P') + V^^^^)ix)] ^^{x), 
where V_^^^^{x) = ^2TTXe^^-f°-f^'r^^^^{x) = ^27tX [j^J^^^^ 



(x) - If 



[x] 



(2.4) 
(2.5) 



and since J'^^f^^^ix) , ^^{x) = 0, we even need not to symmetrize or to pass the operator 

product in the r.h.s. of the HEqs (|2.4|1 into the normal-ordered form, but the careful 
definition of the Heisenberg "currents" V_^(^q,^{x) in Eq. 1)2. 5p is still necessary. On the 
other hand, in order to linearize Eq. p.4|) . we should first determine the dynamics of these 
"currents" . For this purpose we rewrite them in a more convenient abstract form: u = A, 



A 



(2.6) 



and then by means of Eqs. ()2.1() - ()2.6() . we find for these "currents" the formal HEqs 



idQV^^(^<i,){x) = [F_g(vi,)(x) ,ifo(vi,)(t)] + [^-5(*)(a:),i?/(*)(t)] 
= \ ( ^ilAx) \0^E_^{P'^) - S_5(pit)o^j ^_^{x) 



A 



+ (^l^(x) [o^V^^^){x) - y/(^)(x)0^] ^.i:[x)) I (Oi 
which, since oc Of oc /, 7^ ^ ^ (O^) ( [o^ , Of 



A B 



o 



read as idQV_^i^,){x) = \V_^^^){x) ,H^m){t)\ = [yL^(^)(j;) , i?o(*)(i)] • 



(2.7a) 
(2.7b) 

(2.7c) 
(2.7d) 
(2.7e) 



Of course the conditions 1)2. 7d|) are satisfied by any set of the mutually commuting matrices 
O^^O^ . But the final Eq. 1)2. 7ep implies also the absence of possible contributions from 



^Our conventions are: 

,01 _ 



(x , X ), X = t, h — c — 1, dfi, = (9o,9i), metric tensor g 



1, 



1, Dirac conjugate field ^'{(a;) = ^^{x)'y , gamma matrices 7 = ai, 



antisymmetric tensor e ' = ~-e 

7^ — — i(72, 7^ = 7^7^ = 0-3, with (Ti and / as Pauli and unit matrices respectively, and operators 
= -1^1 = -id/dx"-, pit = i ai, E^{P^) = -y^P^ + 7°m^ 
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the Schwinger terms |l23i[ . Only under the assumption, as discussed in the Summary, the 
evolution of operators V_^i^^^{x) is governed by the free one-particle Hamiltonian -ffo(*) (*) 
only leaving just the line (|2.7b|) . A simple main observation of our approach developed 
in ^S], H?], is, that since the contribution responsible for the interaction vanishes, these 
Heisenberg operators (in the Heisenb erg- field representation), as solutions of HEqs (|2.7e|) in 
a weak sense of DM are equivalent to itself in some free-trial-physical-field representation^ 
which, according to the explicit form of (|2.7a|) given below by (|6.1j) . (|6.2|) . we should 
choose as free fields 'tp-^{x) with the same mass from Hamiltonian H2.2|) . In other 
words, from homogeneous HEqs (|2.7ej) with (|2.4j) . (|2.5j) we conclude that the DM T['0(a;)] 
leaves the form of these "currents" invariant up to possible multiplicative constant fixed 
by initial condition at t = equal to unity. That is 

= e^^(*)V_5(^)(x\0)e-*-^(*)* = e^^o(*)V_5(^)(x\0)e-*-^o(*)* < — > (2.8) 



^ e*^owV_5(^)(x\0)e-*^o(^)* = y_5(^)(x), or V-^irm{x) = F_5(^)(x), 

where the last line takes place due to the fact, that just these and only these "currents" 
(|2.5|1 determine the interaction in the HEqs (|2.4|) . so that the latter get a linearized form 

5o^g(x) = -i [E^{P^) + VL^(^)(x)] ^^(x). (2.9) 

Strictly speaking the operator product in the r.h.s. should now be symmetrized or normal- 
ordered. We show below, that the correct result however may be obtained proceeding from 
Eq. (|2.9j) . and argue the only possible weak sense of relation (|2.8|) . 

From a dynamical point of view the suggested linearization condition 1)2. 7e|) . as a free 
equation of motion for Heisenberg "currents" , is a natural and straightforward extension 
of the usual vector current conservation condition. Indeed, supposing for a moment the 
unit matrices in Eq. ()2.6() for all A, = I, we immediately turn Eq. 1)2. 7e|) into the 
well-known divergence identity for the vector current, that is 

i9oJ°^(^)(x) = J°^(^)(x) ,i/o(*)(0 , means 9^, J^^(^)(x) = 0, (2.10) 

oo 

and [Q_^(^)(t),i/(^)(t)] =0, where Q_^^^){t) = j dxV°^(^)(2;\ t). (2.11) 

— oo 

Thus the conservation of the local Heisenberg vector current with corresponding fermionic 
charge Q-^(m){t) is equivalent to the fact, that the time evolution of its zero component 
J^^(\i,){x) is governed by the free part HQ(^q,-^{t) of the total Hamiltonian i7(,j,)(t) only 
(a similar phenomenon in the SG model is described in Ref. [HI). However, since the 
component ^^^^^^(a;) alone does not define entirely the r.h.s. of the equation of motion 



(|2.4j) . then the conservation condition (|2.10|) alone does not enough for the second line in 
Eq. (|2.8j) to exist, whereas the condition (|2.7e|) in common with Eq. (|2.4|) is enough. 

3 Step 2. Functional integral representation of the solution 

The linearized HEqs (|2.9j) allow to write their formal solution in the chosen representation 
of the physical fields 1)2. 8() . as a time-ordered (T-ordered) exponential of different sets of 

^That in fact takes place for all known solvable fermionic models. 
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mutually non-commuting operators, such as = and = —idi, 7^^, 'ijj-^{x) 



^^{x\T)=T 



exp 



*c(x\0), (3.1) 



where the Schrodinger initial conditions at t = are imposed. Following Feynman jl2j . 
the Heisenberg field ((HU) for T > may be reconstructed from the initial field '^^{x^,0) 
using the matrix element of the evolution kernel 3^(^_^), which in turn is expressed as a 
functional integral over the phase space variables x^(r/) and p^{ri) 



00 

^^{x\T)= I dy'y^^_^){T,x'\y\0)^^{y\0), 



yi^_^){T,x^\y\0) = {x^\T 



exp < —I 



dr?(i?5(pi)+y„5(^)(X\7?)) 



(3.2a) 

\y') (3.2b) 



a;l(T)=xl 

J Vx^{ri) j T>p^{r])ex.p li J drj p^ (rj) x^ [r]) 
xi(0)=?/i V 



exp l-i dr] {E^{p^{r])) + {x^{7]),r])) 



(3.2c) 



However, the integrand still contains a "mixture" of different structures, the 7-matrix 
and the field ones 7^^^_^[. . .]. Nevertheless, the explicit form (|2.5|) of the free "currents" 

^-^{tp) (x) = ^27rA 7^ "^-^{tl)) ~ ^ '^-S.ii') prompts the shift of functional variable 
p^{r]) = p^{r]) — ^27rAJ^^^^-| [x^{r]),7]^ to obtain for the matrix element, returning to the 
previous designation p^{r]) l—^ p^{r]) with x'^{r]) = {l,x^{r])), the following expression 



y^^_^^{T,x^\y\0)= I Vx\r])%p_ 

zi(0)=?/i 



exp <-i^27rA jdr]e^ux'^{r])J'^^(^^^ {x^{r]),it]) 



exp —i J dr] {j^p^if]) + ^'^m^ — p^ {ii)x^ {ri)) 



(3.3) 



The most important property of this expression is the splitting of the 7-matrix and the 
field operator orderings, that moreover, are integrated over different functional variables. 
It is a simple matter to see that up to this point the space-time dimension does not play 
any essential role: if the linearization conditions (|2.7d|) . 1)2. 7e|) are satisfied in higher- 
dimensional cases as well the above transformations yield results |16j . similar to Eq. (|3.3j) . 
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4 Step 3. Bosonization and T-ordering 
4.1 Bosonization 

The local vector current of the free massive Dirac field is defined by means of antisym- 
metrization or normal-ordering of a formal field product |3j . Due to the particular property 
of two-dimensional space-time, like 7^7^ = e'^'^j^ and the conservation law similar to (|2.1fl|l 
these free massive currents bosonize to pseudo-scalar fields <l>_g(x) |U, ^\ 



leading to the simple full-derivative representation 

1 d 



where J.^.^Jx) =: V-5(a:)7^V'-c(a;) :, 



IT dr] 



^-5 {x\v),v) 



(4.1) 



(4.2) 



that recasts the first of the time-ordered exponentials in Eq. to the form 



T 



exp I -ig j dr]-^ {x^{rj),r]) 



with g = 2^A/7r. (4.3) 



Thus, the problem of T-ordering of the exponential bilinear in the fermion field is trans- 
formed to the T-ordering of an exponential linear in the pseudo-scalar field $_^(x). This 
is a considerably more simple exercise at least for the free boson case, as is shown in Ap- 
pendix. The necessary relation between the pseudo-scalar fields $_g(x) and free pseudo- 
scalar fields 0_^(x), turning the problem to the free-field case, is found in the next sub- 
section. 



4.2 Reduction to the free massless pseudo-scalar fields 

A simplest way to relate the field <I>_^(x) to the free field (/)_g(x) can be derived from the 
relation between the free massive ip-^{x^,t) and the free massless X-S,{x^^'t) Dirac fields, 
with the formal unitary transformation operator G_^(t) [3], so that for arbitrary t, from 



i-f^df.Tp^^^x) = m_^ip{x), il^d^x-ii{x) = 0, 
one has ^/;_^(x\i) = Gz\{t)x-^i{x^ ,t)G^i^{t), 
as well as J^g(^)(x\t) = G-_\{t)j'^,M{x\t)G_i:{t) 



where f_^i^^){x) =: X-ii{x)-i^X-d^) ■ ■ 

Now we temporarily omit for brevity the index — ^, using 
by means of these Dirac equations the following relations 

[R{t),x{x\t)]=imj\{x\t), for R{t) = {doG{t)} G'^t), 
00 

so, that R{t) = -im j dy^ {xiy\t)x{y\t)) = -R^{t), 



(4.4a) 
(4.4b) 
(4.4c) 
(4.4d) 

m and so on, to obtain 



G{t) = 



exp I J driR{r]) 



G(0), GHt) = G-\t), 



(4.5a) 
(4.5b) 

(4.5c) 
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where the solution ()4.5b|) . (|4.5c|) of Eq. (|4.5aj) in terms of operators xi^^^t) is obtained 
by means of the corresponding equal-time anti-commutation relations of the type (|2.3() . 

The following comments are in order. First of all, we note, that an operator of the 
same type as R{t) is present in the Hamiltonian 1)2. 2|) and should exist under appropriate 
regularization for example in normal-ordered form [HI, |22]. The G(t)-transformation re- 
lates the time evolutions of two different free quantized fields, each of them is defined on 
a corresponding space of states. In fact, it is a Bogoliubov transformation diagonalizing 
the free massive fermionic Hamiltonian in terms of the free massless fermion fields for ar- 
bitrary time t similarly to |4j. The operator G(0) can be used to impose additional initial 
conditions. Because this operator does not change the final results of further calculations, 
we eliminate it further on choosing for simplicity G(0) = /. 



For the free massless fermion current ()4.4d|) the above arguments, originative earlier 
the bosonization rules 1)4. introduce now the free massless pseudo-scalar field (j){x) jl3j 

W(x) (^) = ^ df,(t){x) , x(.3:)x{x) = V cos {(3(t){x)) , (3 = 2^7?, (4.6) 



/vr 

connecting this current and the field by the equation similar to Eq. (|4.2I) . so that both 
the equations amount to the relations between the pseudo-scalar fields only 

<i<^(.-C),.) ,g-. = G-'W*(x',t)GW. (4.7) 

dr] \ ar] I 



where G{t) =7^ 



exp < J dr]R{r]) )■ , R{r]) = -imv / dy'^ cos {f3(p{y^ ,r])), (4.8) 
I 

and the index — ^ is assumed for all quantities. The value of v depends on the renormal- 
ization prescription chosen for the normal-ordering procedures of boson and fermion fields 
[0], jlOj assumed in both sides of Eqs. (|4.6|) . These procedures are however irrelevant for 
our further operator transformations, so we do not specify further the value of v. 

The next comment concerns the properties of the field $(x). These may be easily 
derived substituting into the free field equation for the field (p{x) its formal expression 
from Eq. (|4.7|) . Namely, for x^ = t 

a^9^0(x) = {di - df) (l){x\t) = 0, where ^{x\t) = G{t)^{x\t)G-^{t). (4.9) 

Using here Eqs. (|4.8|) we arrive at the following form of the second-time-derivative 

dicp{x\t) = G{t)[dmx\t)]G~\t) - [doHx\t),R{t)]. (4.10) 

The commutator is easy computed by means of (|4.8() and the canonical commutation 
relation for the free pseudo-scalar field [do(j){x^,t) ,0(y^,t))] = —i6{x^ — y^) as follows 

[5o(/>(x\t),i?(t)] = mi;/? sin (/3(/>(j;\t)) . (4.11) 
Thus the second relation (|4.9|) leads to the SG equation of motion for the field $(x) 

^^,^^'(^[x^,t) = G{t)[d^d^^{x) - mvl3s\ii{P^{x))]G~^(t) = 0, (4.12) 
which is reproduced for (3 = I^/tx by the well-known Lagrangian density [Tj, QZ]; [13 

L^G{x) = \^^^,^{x) d^^{x) - vm [cos (/3$(x)) - 1] . (4.13) 
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4.3 Conversion of the T-exponential into the unordered form 

Now we are ready to get rid of the T-ordering in Eq. (|4..S|) . This problem is simphfied by 
using the well-known operator formula for the T-exponentials ^Hl 



Ua+b{T)=Ua{T)Uc{T), where C{r]) = {UA{r])) B{r])UAiv), 

T 

exp <^ / dr](A{r]) + B{rj) 



(4.14) 



with Ua+b{T)=T 



and so on. 



Identifying here A{ri) = R{r]), B{r]) = —ig(^d(/){x^{r]),ri)/drj), by means of the first Eq. 
(|4.7() and Eq. (|4.8() . we transform the expression on the r.h.s. of Eq. 1)4. 3|) as follows 



exp 



dr] R{r]) - ig 



d(l){x^{rj),r ])' 
dr] 



, with g = 



(4.15) 



According to Eqs. 1)7. 10(1 . ()7.11|) of Appendix, the T-exponential ()4.3() for the case of a 
free (pseudo-) scalar field is independent of the trajectory x^(r/) and depends only on its 
end points 



W{T) = 



exp 



ig / dr] 



d(j){x^ (rj) , rj) 
dr] 



U{T)U-\0). 



where U{T) = exp 



l-{±\T)y exp[-igcPix\T),T)] 



(4.16) 



(4.17) 



that allows to apply again the relation 1)4. 14() . written now from left-to-right, to the last 
7^— exponential in Eq. (|4.15|) . just interchanging the roles of operators only: B(r]) = R(r]), 
A{r]) = —ig (^d(p{x^ {r]) , rj) / dr]) , recasting Eq. 1)4.15(1 into the form 



G~\T)U{T)r^ 



exp 



J dr] U-^{r])R{rj)U{r]) 



u-Ho). 



(4.18) 



From the Eqs. (|4.8|) . (|4.17() . with the locality condition for the field 4>{x), originated in 
Eqs. ()7.H) . (|7.2|) of Appendix, it is a simple matter to see that, when x = x{r]) = (77, x^{r])), 
y = {V, y^), then x{r]) - y = (0, x'^{r]) - y^), 



[cl)ixir])),<Piy)] = -iDoixi7])-y) = 0, U-\r])R{r])U{r]) = R{r]), 



(4.19) 



and the integrand in Eq. ()4.18|) is reduced to R{r]). Therefore, by means of (|4.7() . we 
arrive at the final expression for the T-exponential of Eqs. (|4.3|) = (|4.18|) . as 



G-\T)U{T)G{T)U-\Q) = n{T, 0) exp [-ig^x{T))] exp [i<7$(x(0))] , 



where 0,{T, 0) = exp — i— 



1 



1 



.1(0))' 



(4.20) 
(4.21) 



is a c-number phase factor. It different from unity, for example, if absolute value of the 
velocity 1x^(^)1 during the evolution has crossed the speed of light an odd number of times. 
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5 Step 4. Initial and Heisenberg fields 



Now we are finally able to obtain the solution of the Federbush model [21] • First 
of all we notice that the above obtained expression (|4.2())) for the T-exponential is fully 
independent of the trajectory x^{r]) and factorizes out from the path integral in Eq. H3.3|) 
which is then immediately calculated by definition (|,3.2a|) . (|,S.2c|) as free evolution kernel 
1^'''''* , expressing the time evolution of the another free fermion field with the mass 



y{^-AT, x^\y\0) = n{T, 0) exp [-tg<^^^{x{T))] exp [ig^.^ixiO))] Y^'^^ {T, x^\y\0) 



where Y^°\t, x^\y^ ,0) = 0(r)exp [-iTE^ (P^)] 6{x^ - y^) 



(5.1a) 
(5.1b) 



a;l(T)=a;l 



xi(0)=yi 



exp 



i / dr] (^■j^p^{r]) + "y^m^ — p^{r])x^{r])) 



so, that J dy'Y^''\T,x'\y\0)iJi:{y\0) = iJi:{x\T). 

— oo 

In order to check Eq. ()5.1b|) we need only the following relation 



(5.1c) 



Vx'^{rj)6 { x^ -y^ - / drix^{r]) [. . .] 



Vx\v) [■■■]-- 

xi(0)=j/i 

and the Fourier-representation for the delta-function. Prom the Eqs. (|5.H) above it is a 
simple matter to see, that in order to obtain the correct solution for x^(T) = x^, the initial 
field for x^(0) = y^ should be chosen as follows 

^^{y\0) = exp [- ig^^^{y\0)] ^^{y\ 0), g = 2CXV^, (5.2) 
so, that ^5(x\T) = T[$_^,V'5] = ri(T,0)exp [- ig^^i:{x'^ ,T)] ^^(x\r). (5.3) 

It should be stressed that, as a solution of the linear homogeneous equation ()2.9() . this 
expression may be easily renormalized following Wightman 2_lj, |29j by dividing the r.h.s. 
of Eq. H5.3|) on the vacuum expectation value (0| exp [— ig^_^{x))] |0). The factor 0(r, 0) 
1)4. 21(1 is generated in Eq. (|7.9() of Appendix by non-vanishing contribution of zero modes 
into the dimensionless Pauli- Jordan commutator function (|7.2|1 of the two-dimensional free 
massless (pseudo-) scalar fields. As was shown by Faber et al. [S], ^0], these collective 
modes are not affected by the dynamics and their consecutive removal makes meaningful 
the theory of the free massless two-dimensional (pseudo-) scalar field. With this interpre- 
tation the value |x^(t)| is the constant velocity of the "center of mass system" of these 
zero modes decoupled from everything, and i^{T, 0) = 1 for all T. 



6 Summary 

We have shown, how the solution of the Federbush model may be obtained by direct in- 
tegration of the linearized formal HEqs using the notion of physical fields. These fields 
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are conventionally taking as the Schrodinger ones, because the latter as well as Heisen- 
berg fields form a complete set of fields itself, unlike the asymptotic in (out) fields PH] . 
what is also useful for construction of explicit solutions of the bound-state problem 
j27j . We have revealed the fundamental role of the free massless pseudo-scalar field for 
this integration procedure. We have found, that the solvability of the Federbush model, 
besides the properties of the linearization ()2.8|) and the bosonization H4.1|) . hides in the 
full-derivative expression (|4.2|) for the interaction term with an arbitrary trajectory in the 
time-ordered exponential of the path integral representation for the evolution kernel H3.3|l 
of the Heisenberg field. It is worth to note that the method suggested here may be applied 
to another known solvable two-dimension models like the Thirring and derivative coupling 
models ^Hj) the model in Ref. as well as to some four-dimensional ones |16j . 

Our last comments concern a possible contribution to Eq. (|2.7a|) from the Schwinger 
terms. For the massless Thirring model it may be shown that their invariant contribution 
vanishes independently of the value of the constant in the current commutator. On the 
one hand, it is well-known [H], |2H]) that this constant directly depends on the chosen 
physical field representation - the chosen space of states. On the other hand, it was shown 
[23) . that for the Federbush model its value strongly depends on the chosen point-splitting 
prescription for the Heisenberg vector current in the HEqs ()2.4j) . (|2.5|1 . ()2.7j) . That is why 
in this work we prefer to deal with local currents only. 

Indeed, for these currents, because of 7^ = 7*^7^, the "non-interacting" and "interact- 
ing" parts of the HEqs (|2.7aj) - ()2.7c)) may be transcribed straightforwardly as follows 

idoV_^^^)ix) - [V-^i^)ix),Ho(^^)it)] = 7' <7^/^|i9, {^.^ix)Y^-d^))'^ (6-1 

[y„5(q,) (x) , Hj^^) it)] = 7^ gV^ [ (x) , Hj^^^ (t)] (6.2) 
-IqVt^ J1^(^)(x) 

Matching the independent matrix structures / and 7^ of these two equations we easily 
recognize both vector current and axial current divergence identities if the "interacting" 
part ()6.2)) vanishes identically. Therefore the conservation of the Heisenberg vector current 
1)2. 1U() preserves here the absence of the full contribution (|6.2() from the Schwinger terms 
in the r.h.s. of Eqs. (|2.7aj) - ()2.7c)) . because for the given interaction Hj(^^^{t) ^2.1^ . (|2.2jl 
both commutators in the r.h.s. of Eq. (|6.2() reduce to the one and the same Schwinger 

term in the one and the same commutator J^c/q,) {x) , J^c/^) (y) , and are vanishing 

or not only both together. It means that the bosonization rules like ()4.1() . (|4.6() cannot be 
applied here to the local Heisenberg vector current J^^(^-)(x) directly, unlike 0, [HI, [13] . 

The observed contradiction may be overcomed following Ref. [25 in QED, if we impose 
bosonization rules to the free fermion currents ()4.1I) . (|4.6|) only, and if only weak corre- 
spondence (|2.<S|) between the Heisenberg and free currents take place. Thus the Schwinger 
term does not contribute into the HEqs (|2.7a|) . because it becomes meaningful only af- 
ter the identification transformation ()2.8)) determining the representation space. These 
arguments, besides of all, demonstrate a self-consistence of the suggested linearization 
procedure ((231), ^ ^ UMi with obtained solution JOl) or JO)), ([0[l . 
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Thus we come to the conclusion: the transformation of the "currents" (|2.8)1 . generated 
for the Federbush model by the DM (|5.3|1 with the pseudo-scalar SG field and 
normal-ordered with respect to the Schrodinger free-fermion-physical fields ip±^{x) 

where ^-^{x^,t) 

(see ()4.1(l ) leaving invariant the local form of the both components of vector current 
J-^(^q,)ix) = &s well as conservation law (|2.1U() . changes their algebra, extending 

it by the Schwinger term, whereas for the conserving charges (|2.11|1 

<3-?(*)(0 = Q-m(*> = ^ = (6-5) 

the algebra remains unchanged. 

We are grateful to Prof. Vail A.N. and Prof. Leviant V.M. for useful discussions, and 
to the referee of JNMP for suggesting improvements of the manuscript. 



exp — 

( oo \ 



(6.3) 
(6.4) 



7 Appendix 

Let us consider more general time-ordered exponential of the total derivative of the free 
massive (pseudo) scalar field (p{r]) = gcp{x^{r]),r]), taken on an arbitrary trajectory x^(ry), 
where its commutator function Am{t\r]) is expressed via the Bessel function ^o{z) as 



[ip{t),ip{r])] = A^{t\rj) = -ig^D^ {x{t) - x{r])) 



Dmix) 



27ri 



(2vr)2 



-i(kx) 



e{xo) 



e (x2) Jo 



mV X 



(7.1) 
(7.2) 



Since this commutator is a c-number, by repeated application of the well known formula: 
exp(74) exp(i?) = exp(^ -|- B) exp({l/2}[^, B]), one can obtain for the T— exponential 



exp < 



dt 



dtf{t) 



dt 



■ exp 





T V 



exp{-iip{T)} exp {MO)} exp { - -A™(r|0) 



(7.3) 



and transcribes it, by means of the locality condition (|4.19|) and the next equality 



r rj '1' 

j drj j dr ||^^A„(r?|r)| = A„(T|r) - A„(r|0) - j 



drj 



d 

d^ 



Am(?7|r] 



T=r] 
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as 



exp {— i99(T)} exp {i(^(0)} exp 



dri 



d_ 
dr 



Am(??|r) 



T=T\ 



(7.4) 



The function Am,(^|T) is discontinuous on the hght-cone, at r] = t, and its derivative has 
a singularity there. Therefore the expression for the integral in the last exponential is 
indefinite and should be redefined as the following limit with a +0 of 



dr] 



dr 



T=ri—a 



^ I dr] 27ri 



(2^ 



e(A;°)5(A:2 - m^) 



(7.5) 



-i{k-{x{ri)-x{T))) 



T=rj — CF 



For the latter derivative, using the Taylor expansion of the trajectory x^(r/ — a) with 
respect to o" > 0, and changing the momentum variables as = ak^, we obtain 





' d 


1 d ' 


T=r)—a 


dr] 





dr 

After these changes the integral (|7.5jl becomes 

T 



(7.6) 



7=1 
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dr] 





I. 
2i 



dr] 



d Id 

dr] a 

d 1 d 

dr] a 



2-Ki 



(2^) 
e(7)0 (7' 



2^2x -i7(g"-g'AM'?)) 



(7.7) 



7=1 



l-[± 



' iv)) '] ) Jo iam^^l-{x\ri)y 



7=1 



Because Jq{z) = 1 + O(z^), it is easy to see, that the mass dependance of (j7.2|) leads here 
just to terms of order a. So, as cr ^ +0, the direct substitution 5{q^ — a'^m?) ^ 5{q^), 
gives for the integral (|7.7jl 



d_ _ ld_ 

dr] a 57 



(7.8) 



7=1 



The integrand in ()7.8|) does not depend on 7 at all near the point 7 = 1, and its cor- 
responding derivative vanishes, eliminating the contribution singular at cr ^ +0. The 
remaining integral on drj is taken immediately, leading to the expression for the limit 



1 

lim / 



dr] 



d_ 
d^ 



Am(r/|r) 



T=r)—a 



2i 



{e[l-{x\T))'] -0[l-{i\of]]. (7.9) 
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As a result, the time-ordered exponential (|7.3|) takes the following form 




(7.10) 



(7.11) 



This result is of fundamental importance for all computations, because it depends on 
the values of an unknown trajectory x^{7]) and velocity x^{r]) only at the initial and final 
points, t = and t = T. Moreover, the phase factor from (|7.4I) . 1)7.91) in 1)7. 1U() becomes the 
same for both the massive and massless cases, that, on the one hand, corresponds to mass 
independence of the light-cone behavior, but on the other hand, due to the dimension- 
less nature of the commutator function (|7.2() . represents by the construction H7.5|) ~ H7.8|) 
the non-trivial infrared properties of the free massless (pseudo-) scalar field, such as the 
existence and the non- vanishing contribution of its zero modes [0] , ^01 • 
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